Universal properties of Dark Matter halos 
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We discuss the universal relation between density and size of observed Dark Matter halos that was recently 
shown to hold on a wide range of scales, from dwarf galaxies to galaxy clusters. Predictions of ACDM N-body 
simulations are consistent with this relation. We demonstrate that this property of ACDM can be understood 
analytically in the secondary infall model. Qualitative understanding given by this model provides a new way to 
predict which deviations from ACDM or large-scale modifications of gravity can affect universal behavior and, 
therefore, to constrain them observationally. 

PACS numbers: 95.35.+d, 98.62.Gq 
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Introduction. The nature of Dark Matter (DM) is one of 
the most fundamental puzzles in modern physics (see e.g. [1] 
for a recent review). Resolution of this puzzle may poten- 
tially shed a light on many fundamental physical issues, in- 
cluding particle physics beyond the Standard Model, large- 
scale properties of gravity, etc. The simplest Cold Dark Mat- 
ter (CDM) model is rather successful in explaining cosmolog- 
ical and astronomical observational data at large distances |2]. 
More complicated models, discussed in the literature, include 
"warm" DM models [3], or models with modifications of 
large-scale properties of gravity [4] or even of Newtonian dy- 
namics [5]. Search for deviations from the CDM predictions 
provides an important tool to constrain the properties of un- 
known DM particles and to test fundamental laws of gravity. 

Comparison of predictions of DM models with observa- 
tional data is complicated by the fact that information about 
the properties of DM halos is usually derived from the study 
of dynamics of "luminous tracers": stars and gas particles (the 
only exception being gravitational lensing method). Physics 
of the luminous baryonic matter in galaxies and galaxy clus- 
ters is very nontrivial to be fully accounted for in theoretical 
and numerical modeling. Thus, it is crucial to find the prop- 
erties of DM distributions that do not depend on the physics 
of the baryonic matter and, therefore, could be directly com- 
pared with predictions of pure DM models. Below, we argue 
that recently discovered relation between characteristic den- 
sity and size of the DM halos of galaxies and galaxy clusters 
may provide an example of such a "universal" relation, and 
suggest its analytical understanding within the CDM model. 

Universal properties of DM halos. It was discussed for 
a long time (see |6r|8J] and references therein) that sizes and 
densities of observed galactic DM halos obey a universal scal- 
ing relation. Namely, the DM central density pc, averaged 
over a region with the characteristic size tq of the halo, scales 
roughly as r^ 1 . This scaling law might appear surprising, if 
one takes into account quite different histories of formation of 
various galaxies and their different types. No clear qualitative 
understanding of the origin of "dilution" of the density of DM 
halos with the increase of the halo size was proposed so far. 



A different scaling relation between parameters of DM dis- 
tribution was suggested in |01 for still larger range of distance 
and mass scales. The authors of |0] collected published data 
on the DM distributions in 289 objects ranging from dwarf 
spheroidal galaxies (dSphs) to galaxy clusters. The obser- 
vational data were fit by different density profiles: Navarro- 
Frenk- White (NFW) [ 10], (pseudo)-isothermal (ISO), Burkert 
(BURK) 111 111 and on average 3 density profiles were collected 
for each object. It was found that the average "column den- 
sity" of DM halos 
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slowly grows with the increase of the mass of the halo as 
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where M200 1S the measure of the halo mass (it is the mass 
contained within a sphere of the radius -R200 where the aver- 
age DM density equals to 200p cr i t with p crlt being the critical 
density of the Universe today). The parameter rc in (Q~|i is 
taken to be r s of the NFW profile that fits the observed DM 
distribution. If the same data (e.g. rotation curves in galaxies 
or temperature profiles of galaxy clusters) were fitted using a 
different DM profile, the scale rc may be uniquely expressed 
in terms of the parameters of that profile. The column den- 
sity (Q~|i is insensitive to the type of DM density profile used to 
fit the same observational data (see Q for details). 

FromEq. (Q~|i it follows that S oc pc r c- Therefore, for cored 
profiles (ISO or BURK) the relation (O is a generalization of 
the previously considered relation pc r c — const 00]. How- 
ever, the quantity S is more universal, as it is defined for any 
(not necessarily cored) DM profile. Therefore, in 19J] S was 
derived for any type of DM density profile, used by observers 
to describe the DM distribution in halos in the broad range 
of sizes (0.2 kpc < r c < 2.5 Mpc) and masses (1O 8 M < 
Af 2 oo < 10 16 M©). The relation p c r c ~ 141M© pc" 2 , de- 
rived in |01 for galaxies, and S — 2OOM0 pc~ 2 , proposed in 
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the Ref. lfl2ll for galaxies and galaxy clusters, could serve as 
approximations of relation (f2j) at low- and high-mass ends of 
the data. This is clear from Fig. Q]in which the data from [9] 
are presented together with the fit © (solid line). 

The relation (ffj) is derived from the analysis of data on 
DM halos of structures with significantly different physics of 
the luminous matter, (from stars in dwarf galaxies through 
gaseous disks at the outskirts of spiral galaxies to hot intra- 
cluster medium gas in galaxy clusters). The fact that the prop- 
erties of DM distributions in all these structures follow the 
same scaling relation allows to conjecture that the observed 
scaling is a "baryon-independent" characteristics of DM ha- 
los. 

A strong argument in favor of such a conjecture is given by 
agreement of the relation (O derived from observational data 
with a similar relation derived from pure DM N-body simula- 
tions, performed in the framework of the ACDM model 111 311 . 
The distributions of DM particles in simulated halos are fitted 
by a density profile, from which S is computed using Eq. (Q]l. 
As discussed in jjj, not only the trend S ~ M ' 2 , but also 
the scatter of individual objects around the average (the pink 
shaded region in the Fig. Q] represents 3c scatter of the data 
from the simulations of Ref. [13]) and a different slope for 
dwarf satellite galaxies (gray dashed line) are reproduced ex- 
tremely well. The quantity S is computed using the best-fit 
models obtained by different groups of observers for each ob- 
ject and compared directly with the results of numerical sim- 
ulations, without questioning the suggested values of M and 
S. The agreement with N-body simulations confirms that, de- 
spite different nature and quality of the observational data, the 
relation|2]is a physical effect rather than a consequence of an 
observational bias. Our goal here is to explain this property of 
the ACDM model. 

Secondary infall model. In spite of apparent simplic- 
ity of relation (fJJ no qualitative explanation and/or analytical 
derivation of the observed decrease of the average halo density 
with the increase of the halo mass or size has been proposed so 
far. In what follows we show that this effect could be readily 
understood within a semianalytical model of structure forma- 
tion, known as "secondary infall" model. Within this model, 
the column density S is expected to change slowly with the 
increase of the overall halo mass, S ~ M£ alo , with n < 1/3. 
Under some simplifying assumptions (see below), the upper 
limit k — 1 /3 is achieved in structures which continue to grow 
at present. 

The secondary infall model 1 15 - 2(J (see also 1 21 ] for a re- 



cent review) provides an analytic description of the growth 
of the mass M(t) and size R(t) of DM halos in the course 
of spherically symmetric accretion from cosmological matter 
flow. One assumes that initial overdensity at the time ti is dis- 
tributed in a spherically symmetric way so that the initial mass 
distribution depends only on the distance r from the inhomo- 
geneity center 
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FIG. 1 : Column density S as a function of halo mass M200 • The gray 
solid line is the relation (fj}. coinciding with predictions of N-body 
simulations fl3ll . using the WMAP5 |2] cosmological parameters. 
The shaded region shows the 3<j scatter in the simulation data. The 
vertical lines indicate the mass range probed by simulations. The 
gray dashed line shows the results of the Aquarius simulation fl4ll 
for sub-haloes. The black dashed-dotted line follows from the infall 
model. 



where pm(U) = (fwrGfl) 1 is the cosmological matter den- 
sity and 8m(r) is distance-dependent mass excess. 1 

The entire DM distribution is split into thin shells parame- 
terized by their initial radii = r(tj). Evolution of the radius 
r(t) of each shell is governed by the Newtonian dynamics 
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where U(r) is the gravitational potential and m(t) is the mass 
inside the radius r(i). Initial velocities of the shells are as- 
sumed to follow Hubble law r(ti) = H{ti)ri, with H(t) be- 
ing the expansion rate of the Universe. 

If <Jm, (r,) > expansion of the shell with given n slows 
down faster than it would do in the absence of central overden- 
sity. A straightforward calculation allows to find the time t* at 
which r(i«) = and the shell reaches maximum (turnaround) 
radius i?* = r(<*) 



n 2 Rl 
8Gm, 



1/2 



' 5mi(ri) 



(5) 



At any given moment of time t one can define the boundary 
of the spherical halo R(t) as the turnaround radius of a shell 
for which t = i*. From Eq.© the mass within this radius, 
M(t) = mi, is related to R(t) via 



R(t) = 
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Continuous growth of the size and mass of the spherical ha- 
los is similar to spherically symmetric accretion onto a central 
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We neglect A-term for sufficiently early t . 
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mass concentration. In the accretion theory, the "sphere of in- 
fluence" of a body with the mass M is defined as the distance 
range within which the gravitational potential U = GM/R 
is larger than the specific kinetic energies of accreting par- 
ticles. In cosmological settings, the given kinetic energy of 
particles is that of the Hubble flow, K = (H(t)R) 2 /2. For 
each R it decreases with time. This leads to the continuous 
increase of the sphere of influence of the halo with the mass 
M: R oc M 1 / 3 /^ 2 / 3 and, respectively, to the growth of the 
halo mass M, as Eq. (O demonstrates. The accretion from the 
cosmological flow could stop if the "sphere of influence" of a 
given DM halo starts to intersect with the "spheres of influ- 
ence" of the neighboring halos. 

At late stages of evolution of the Universe, the rate of ac- 
cretion onto a DM halo could diminish due to the influence of 
the cosmological constant onto the dynamics of cosmological 
expansion in the direct vicinity of the halo. Modification of 
the gravitational potential entering Eq. 
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by a cosmological constant (with the energy density p\) leads 
to a modification of the relation (f5]l between the i?* and to 



A = Anp A R 3 J(3M) (8) 



where /(A) = (2/tt) dx/y/{ X - 1 - 1) + \{x 2 - 1). For 
the fixed t t , the relation ([8]l has the form A 1 / 2 /(A) = const. 
This means that, for the fixed f„, A = const and, therefore, 
the relation between R* and M, implied by ©, still holds, 
although with « 20% different normalization. 

Once the shell passes the turnaround radius, it starts to os- 
cillate around the center of the halo. The shell will cross other 
shells and its dynamics will become dependent on their move- 
ment. The amplitude of oscillations decreases with time, as 
long as accretion from the cosmological matter flow contin- 
ues. The average (over the shells) radial density profile of the 
halo changes with time. In general, no analytic solution could 
be found for the dynamics of each shell and the average halo 
density profile. However, the study of the halo growth sim- 
plifies if one considers initial perturbations with power-law 
profiles 
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where e is the power-law index and M is a normalization 
constant (see B19L Sec HI.D] for the discussion of the relation 
of the value of e and initial spectrum of density perturbations). 
In this case the halo density profie evolves in a "self-similar" 
manner: its shape can be expressed as 
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The function F(x) is the "self-similar" (time and mass scale 
independent) DM radial density profile with a characteristic 



radius xc- Self-similar solutions exist also in a more general 
infall models, considered in lfl9Tl . One can introduce for each 
shell an angular momentum £, while still keeping the problem 
spherically symmetric fl^l . If the angular momentum I = 
ji? 2 /i* with j = const, the initial conditions (0 still lead 
to a self-similar solution. The detailed shape of the function 
F(x) can be found numerically. Its asymptotic behavior is 
given by H3] 
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The property of self-similarity allows to derive the behavior 
of the column density S oc pc^c = pc%cR{t) as a function 
of M: 
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where in the second relation we use the Eq. (O and xq = 

const. 

Equation (fT2l expresses the universal scaling relation be- 
tween halo column density and mass. It is a direct analog of 
relation (O within the secondary infall model. Notice, that 
it provides an explanation, at least at the qualitative level, of 
the observed slow growth of the DM column density with the 
increase of the halo mass. At the same time, the Eq. (fT2l con- 
tains, apart from the mass M, an additional parameter t, the 
time during which the halo accretes from the cosmological 
matter flow. 

Notice that Eq. (0 (or its analog ([8])) implies that all struc- 
tures which accrete up to the time t have the same average 
density p R (t) = M(t)/ [(4tt / '3) i? 3 (*)] equal to const x 
PM(t) (the cosmological density at the time t). For cosmo- 
logically close objects (including those, analyzed in 19[]) ac- 
cretion still continues, i.e. t = to (the present age of the 
Universe). If accretion onto a halo has terminated at time 
ifinal < (this is the case, for example, for the dwarf 
spheroidal satellites of the Milky Way), then according to re- 
lation (Tl2l) the DM column density of such objects will be 
higher for the same mass. 

For the density profile of the form ( fTOl ITTb the DM column 
density (Q]) is given by 
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where a is a numerical coefficient which could be found 
from numerically calculated profiles F(x). Choosing a typ- 
ical value for xc ~ 0.02 [19], one obtains S = 133M©/pc 2 
for M ~ 10 12 Af o . The relation O (for M = M 200 ) with 
this normalization is shown by the dashed-dotted line in Fig.Q] 
Discussion. The universal scaling relation between the 
mass and column density of DM halos (Tl2l . derived within 
the simple semianalytical model of self-similar secondary in- 
fall provides a reasonable description of the data presented in 
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Fig. Q] The scaling S ~ M 1 / 3 is determined mostly by two 
properties of the model: ( i) Average density inside the sphere 
of the radius R(t) is the same for all DM halos, independently 
of their mass. It is equal, up to a constant, to the cosmological 
matter density at the moment t, see Eq. © and ([8]). (ii) Radial 
density profiles of all halos are approximately the same (the 
self-similarity condition (fTUli). independently of their mass. 

The statement (i) is just a general consequence of the 
physics of accretion from cosmological matter flow. The 
self-similarity property (ii) is a simplifying assumption. Al- 
though it is present in the simplest solutions of the infall model 
(see 11910 and the data suggest that the deviation from the 
self-similarity are not so big, it may or may not hold in re- 
alistic structure formation. Validity of this assumption could 
be verified via comparison with N-body simulations of struc- 
ture formation. In the simulations, density profiles of halos are 
usually fitted by two-parametric NFW profiles, which may be 
characterized, for example, by the scales re = r s and i?200- 
The NFW profiles which arise in the simulations are known 
(see e.g. B13IP to satisfy the relation 

C200 oc M 2 - °-\ (14) 

where C200 = R20o/ r s oc 1/xc- The dependence of ip upon 
mass scale is rather weak and the density profiles predicted 
by N-body simulations do form, to a good approximation, a 
self-similar subset of all NFW profiles. 

The small deviation from self-similarity, described by the 
Equation ( [Pil l, explains a somewhat better fit to the data 
provided by the N-body simulations. Expressing xc from 
Eq. (0 and substituting it into ([T3]l one finds S ~ M 23 , 
rather than S ~ M 1 / 3 , which is closer to the S{M) scaling 
relation © derived from the observational data. 

In the frameworks of the infall models deviations from 
self-similarity arise from slight dependence of e in initial 
conditions (0 and/or the angular momentum I on the mass 
scale B19I1 . Finally, dependence of the S(M) scaling (TT3T > on 
tfinai introduces additional correction for structures in which 
iflnal < to- More precise understanding of the deviations 
from the simplest approximation, discussed above, together 
with better quality of the data, would allow to observationally 
constrain the details of cosmological DM halo formation pro- 
cess. The detailed analysis of the small deviations from the 
self-similar infall model is, however, beyond the scope of the 
present Letter. 

Conclusions. The universal S—M scaling in the DM halos, 
found in the observational data and in pure DM N-body sim- 
ulations can be analytically understood in the secondary infall 
model. It seems to be insensitive to the presence of baryons 
and to the details of DM density distributions. This shows that 
this relation has pure DM origin (in contrast with its interpre- 
tation in favor of Modified Newtonian Dynamics, as discussed 
e.g. in H). 

Qualitative understanding of the scaling between the den- 



sity and size of DM halos discussed in this Letter, opens a pos- 
sibility to estimate expected modification of the scaling over 
a broad parameter space of alternative models of DM (such as 
e-g- OD) and/or large-scale modifications of gravity (see 12211 ') 
and, therefore, distinguish between different classes of such 
models. In this respect, the simplified analytical approach of 
the secondary infall model provides a valuable alternative to 
resource consuming N-body simulations. 

The predictions of the secondary infall model can be tested 
in the Local Group |2(ill . using e.g. the data on the local Hub- 
ble flow l23ll . Our results also motivate dedicated astronomi- 
cal search, especially for the isolated galaxies with the masses 
below 10 10 Mq (two pink downward triangles in Fig. [T} and 
systematic studies of DM content in galaxy clusters, as well 
as a uniform analysis of all observations collected in |9j]. This 
would allow to reduce systematic errors in the data, determine 
more precisely the slope in relation (fJJ and find possible devi- 
ations from the pure CDM model. 
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